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ABSTRACT 

MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is an integral-field spectroscopic 
survey that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). 
MaNGA’s 17 pluggable optical fiber-bundle integral field units (IFUs) will observe a sample of 10,000 
nearby galaxies distributed throughout the SDSS imaging footprint (focusing particularly on the North 
Galactic Cap). In each pointing these IFUs are deployed across a 3° field; they yield spectral coverage 
3600-10,300 A at a typical resolution R ~ 2000, and sample the sky with 2” diameter fiber apertures 
with a total bundle fill factor of 56%. Observing over such a large field and range of wavelengths 
is particularly challenging for obtaining uniform and integral spatial coverage and resolution at all 
wavelengths and across each entire fiber array. Data quality is affected by the IFU construction 
technique, chromatic and field differential refraction, the adopted dithering strategy, and many other 
effects. We use numerical simulations to constrain the hardware design and observing strategy for the 
survey with the aim of ensuring consistent data quality that meets the survey science requirements 
while permitting maximum observational flexibility. We find that MaNGA science goals are best 
achieved with IFUs composed of a regular hexagonal grid of optical fibers with rms displacement of 
5 /mi or less from their nominal packing position; this goal is met by the MaNGA hardware, which 
achieves 3 /mi rms fiber placement. We further show that MaNGA observations are best obtained in 
sets of three 15-minute exposures dithered along the vertices of a 1.44 arcsec equilateral triangle; these 
sets form the minimum observational unit, and are repeated as needed to achieve a combined signal- 
to-noise ratio of 5 A -1 per fiber in the r-band continuum at a surface brightness of 23 AB arcsec -2 . 

In order to ensure uniform coverage and delivered image quality, we require that the exposures in a 
given set be obtained within a 60 minute interval of each other in hour angle, and that all exposures 
be obtained at airmass <1.2 (i.e., within 1-3 hours of transit depending on the declination of a given 
field). 

Subject headings: atmospheric effects, methods: observational, surveys: galaxies, techniques: imaging 
spectroscopy 
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1. INTRODUCTION 

Integral field spectroscopy (IFS) at optical and infrared 
wavelengths is among the most significant developments 
in modern observations of galaxies at all redshifts be¬ 
cause it combines the benefits of two-dimensional pho¬ 
tometric analysis with physical diagnostics of baryon 


composition and kinematics (e.g.. lEmsellem et al 
Law et ah] [2009 ; Bcrs hadv all[2010b iSa n chcz et al 


Wcijmans et al 


2004; 


2012 


2014 iFabricius et al.1 l201^b Recent 


advances now enable m ulti-object IFS with instru- 
ments such as SA MI (|Croom et all 120121). KMO S 
dSharples et al.l[20T3h . and MaNGA ( Drorv et al .1120loll . 
As a part of the 4th generation of the Sloan Digital 
Sky Survey (SDSS-IV), the MaNGA (Mappin g Nearby 
Galaxies at APO) project (jBundv et all 12015 1 bundles 
fibers from the BOSS ( Baryon Osci llation Spectroscopic 
Survey) spectrograph (]Smee et al.l I2013T ) into integral- 
field units to obtain spatially resolved optical spec¬ 
troscopy of 10,000 nearby galaxies over a 6 year survey. 
Early results obtained with prototype MaNGA hard¬ 
ware (jBelfiore et alJl2015t iLi et' al.ll2015l : iWilkinson et al.1 
l2015fj demonstrate the richness of the data for exploring 
the stellar and gas composition. 

Because current large-format detectors lack energy res¬ 
olution throughout most of the electromagnetic spec¬ 
trum, IFS has adopted a range of technical approaches 
to down-selecting and formatting a subset of the three- 
dimensional data cube of wavelength and spatial po¬ 
sition onto a two-dimensional detector array. These 
approaches yield different, science-driven trades in the 
data-cube sampling. Simultaneous and integral cover¬ 
age of the spatial field is desirable and achieved by a 
number of instruments us i ng lenslets (e.g., S AURON, 
OSIRIS; IBacon et, al.ll200lt lLarkin et al.ll2003D or image 
sheers (e.g., SINF ONI, MUSE: lEisenhauer et al.l 120031 
IBacon et al.ll2010h . However, the two current wide-field, 
multi-object, IFS instruments-SAMI and MaNGA-use 
bare-fiber arrays to minimize cost while maximizing flex¬ 
ibility and patrol area, but at the penalty of not achiev¬ 
ing truly integral spatial coverage at any one time. This 
shortfall can be overcome by careful attention to the in¬ 
terplay of the hardware design of the fiber bundles and 
the observing strategy. 

The most immediate challenge is that the MaNGA 
fiber bundle, composed of circular apertures with large 
interstitial gaps that significantly undersample the PSF 
at the focal plane of the telescope, has a non-uniform re¬ 
sponse across each IFU. This means that (under most 
techniques for the reconstruction of images from the 
data) the appearance of objects that are small with re¬ 
spect to the fiber size (e.g., AGN or Hu regions) can 
vary across an IFU. The reconstructed image of such un¬ 
resolved objects can either look small and circular (if the 
object was centered on a single fiber), large and circular 
(if the object was centered in the interstitial gap between 
three fibers), highly elongated (if the object was cen¬ 
tered midway between two fibers), along with any range 
of shapes in between. 

This is highly undesirable from a science standpoint, 
and therefore typ i cal fiber-bundle IF U surveys (e.g., 
iCroom et al.ll2014 (Sanchez et al.l 120121) dither their ob¬ 
servations. Small dithers of a fraction of the fiber spac¬ 
ing sample the missing points in the image plane and 


allow reconstructed images based on multiple, dithered 
exposures to achieve fairly uniform and integral spatial 
coverage. 

This dithering is complicated by atmospheric refrac¬ 
tion however, especially given the extremely wide spa¬ 
tial and spectral coverage of MaNGA. Chromatic dif¬ 
ferential refraction over the MaNGA wavelength range 
(AA 3600 — 10300 A) can be comparable to the diameter 
of individual fibers, and field differential refraction (from 
variation in the amount and direction of refraction over 
the 3° field of an SDSS plugplate) contributes similarly. 
These effects combine to degrade the effectiveness of a 
regular dithering scheme in sampling the image plane. 

This paper presents simulations that explore the im¬ 
pact of these effects on the expected MaNGA data qual¬ 
ity, and thereby constrain the hardware design and ob¬ 
serving strategy for the survey. In )J2]we give an overview 
of the SDSS 2.5m telescope and plugplate system, along 
with a brief description of the MaNGA legacy hardware 
and IFU ferrule designs considered for the survey. We 
describe the basic design considerations for the survey in 
521 Using the science requirements summarized in 1 )3.11 
typical integration time s set by the read noise character¬ 
istics of our detectors (1 )3.21) . and numerical simulations 
( D3.3I) we motivate the need for dithered observations and 
regular hexagonal packing of the IFU fiber bundles, cul¬ 
minating in a baseline hardware design and observing 
strategy described in 1 )3.41 This baseline observing strat¬ 
egy is significantly complicated by atmospheric differen¬ 
tial refraction, and we discuss the impact of chromatic 
and field differential refraction on our data quality in 1)4.11 
and 14.21 respectively, defining a uniformity statistic D to 
describe the data quality in lj5l Using the D statistic we 
formulate our final observing strategy in terms of visi¬ 
bility windows in 1JH1 noting a few additional practical 
considerations (e.g., dithering accuracy and IFU bundle 
rotation) in )J7] We summarize our conclusions in lj8] 

2. OBSERVATORY AND HARDWARE OVERVIEW 
2.1. Observatory and Legacy Hardware 

MaNGA opera tes on the SDSS 2.5-m telescope 
(jGunn et al.ll2006l) located at Apache Point Observatory 
(APO; latitude cf) = +32° 46'49"). The telescope is a 
modified Ritchey-Chretien with alt-az mount that is de¬ 
signed with an interchangeable cartridge system that can 
be installed at the Cassegrain focus. T he MaNGA hard - 
ware is described in greater detail bv lDrorv et al.l ( 2015 1: 
here we briefly review the major salient features of the 
system. 

MaNGA has 6 cartridges, each of which contains a 
plugplate with a field of view ~ 3° in diameter that has 
been pre-drilled with holes corresponding to the locations 
of target galaxies into which optical fibers and IFUs can 
be plugged each day in preparation for a night of observ¬ 
ing. These plates are fixed at zero degrees position angle 
(i.e., the on-sky orientation of the telescope focal plane 
coordinate reference frame is fixed). 

Each MaNGA cartridge has a total of 1423 fibers (709 
on spectrograph 1, 714 on spectrograph 2), correspond¬ 
ing to 17 science IFUs ranging in size from 19 to 127 
fibers (12.5 — 32.5 arcsec diameter; 1247 fibers total), 
twelve 7-fiber mini-bundles used for spectrophotometic 
calibration (84 fibers total; see Yan et al. in prep), and 
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92 single fibers used for sky subtraction that can be de¬ 
ployed within a 14’ radius of their associated IFU har¬ 
ness 03 Each IFU has its rotation fixed using alignment 
pins in the ferrules that plug into corresponding align¬ 
ment holes located a short distance West of each target 
galaxy. _ 

These optical fibers feed the twin BOSS dDawson et al.l 
120131 Bary onic Oscillation S pectroscopic Survey,) spec¬ 
trographs dSmee et al.ll2013ll . The collimated beams in 
each spectrograph are split with a dichroic and feed a 
blue (AA3600 — 6000 A) and red camera (AA6000 — 10300 
A). The blue cameras use blue-sensitive 4k x 4k e2V 
CCDs while the red cameras use 4k x 4k fully-depleted 
LBNL CCDs; all cameras have 15/im pixels. Spectral res¬ 
olution varies with wavelength from R = X/6X ~ 1400 
at 3600 A to i? ~ 2000 at 6000 A (blue channel), and 
R ~ 1800 at 6000 A to R ~ 22 00 at 10300 A (red 
channel; see Fig. 36 of lSmee et al.1120131 ). Spectra from 
each of these four cameras are extracted and processed 
through sky subtraction, spectrophotometric calibration, 
astrometric registration, and reconstructed into three- 
dimensional data cubes using a software pipeline (Law 
et al. in prep) descen ded from tha t prev iously used for 
BOSS (idlspec2d; see lBolton et al.l T2012L Schlegel et al. 
in prep) 

The telescope guider system is optimized for a wave¬ 
length of ~ 5500 A and uses endoscopic fibers inserted 
into 16 holes in each plugplate corresponding to the lo¬ 
cations of bright guide stars. These endoscopic fibers 
produce images of the guide stars on a guider camera, 
and the guider actively adjusts the focus, scale, rota¬ 
tion, and offset of the telescope focal plane to track these 
stars through varying weather conditions and observing 
angles. 


2.2. IFU Ferrule Design 

The ability of an IFU fiber bundle to deliver good, re¬ 
peatable, and uniform image quality depends most fun¬ 
damentally on the arrangement of fibers within the bun¬ 
dle; while dithering 1 113.3.211 . differential refraction ($4|), 
and other considerations are important, the fiber place¬ 
ment sets the basis for the sampling regularity of the 
entire survey. 

As described bv lDrorv et al.l (12015D . the MaNGA fibers 
have an inner light-sensitive core diameter (ID) of 120 /im 
(corresponding to 2.0 arcsec in the telescope focal plane) 
and an outer diameter (OD) of 151.0 ±0.5 fin i with their 
protective buffers and cladding. We originally considered 
two kinds of fiber bundles for MaNGA, as illustrated in 
Figure |T| The first was a circular bundle of fibers that 
maximizes the filling factor of light-sensitive fiber cores 
relative to the total IFU footprint by chemically stripping 
the protective buffers from the en ds of each fiber. As de¬ 
velope d for the SAMI survey by iBland-Hawthorn et all 
(1201011 . these ‘Sydney-style’ bundles maximize the effec¬ 
tive filling factor at the cost of decreased fiber throughput 
due to focal ratio degradation (FRD), greater fragility of 
the glass cores, and irregular fiber packing due to the 

14 The physical size of the hardware components also defines 
a minimum-distance exclusion zone around each plugged object. 
These exclusion distances are 116” (7 mm), 89” (5.35 mm), and 
62” (3.7 mm) for IFU-IFU, IFU-sky, and sky-sky fiber placement 
respectively. 
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Fig. 1.— Fiber bundle designs considered for MaNGA (white 
regions represent live fiber cores). The left-hand panel shows a 
127-fiber bundle for which the fibers are arranged in a regular 
hexagonal array (i.e., the final MaNGA IFU design; shown here 
is as-built harness ma024); the right-hand panel shows an exam¬ 
ple bundle of 61 fibers in a circular packing arrangement based on 
that adopted by the S AMI team for use at the Australian Astro¬ 
nomical Observatory dCroom et all [20121 , compare their Fig. 3). 
Although the hexagonal arrangement of fibers has greater regu¬ 
larity, the circular arrangement has greater effective filling factor 
since the protective buffers are stripped. 


circular ferrule geometry. Based on the numerical per- 
formance simulations described in H3.3.3I we prototyped 
(and ultimately chose to adopt) a second style of fiber 
bundle composed of a regular arrangement of buffered 
fibers within a tapered hexagonal ferrule for which we 
pioneered a novel construction technique (see details in 
iDrorv et al.ll2015 ). While reaching lower effective filling 
factor, this technique improves fiber throughput JI3 de¬ 
creases breakageFu and (by virtue of its hexagonal geom¬ 
etry) permits extremely regular fiber placement within 
each IFU. 

The theoretical effective fiber packing density of the 
hexagonal IFUs can be defined as the ratio of the total 
fiber core area (A cor e) to the area of the hexagon circum¬ 
scribing the fiber bundle (A^ex); where: 


Ahex = -^-d 2 (v / 3VR ± l) 2 


( 1 ) 


Acore — 7T ^^ (1 ± 3Ar(Ar ± 1)) (2) 

Here d = 151 fim is the outer diameter of an individual 
fiber, t = 15.5 fim is the thickness of the fiber buffer and 
cladding, and Ar is the number of ‘rings’ in the bundle 
(Ar = 2 for a 19-fiber IFU, and Ar = 6 for a 127- 
fiber IFU). In Figure [2] we plot the effective filling factor 
/ = Acore/Ahex as a function of the buffer thickness t. 
In accord with these predictions, the prototype circular 

15 A conservative e stimate can be made by comparin g Figure 4 
of ICroom et al.l d20121 '> to Figure 11 of IDrorv et al.l (120151 '): MaNGA 
achieves 95 ±1 % throughput with an exit f-ratio of f/4 for fibers 
fed at f/5. In contrast, the original SAMI bundles achieved 50- 
75% throughput with an exit f-ratio of f/3.15 fed at f/3.4. We note 
that the FRD of even t he second-generation SAMI bundles (Fig. 
5 of IBrvant ct alj 12014 1 is sufficiently large that it would require 
our optics to be 40% larger in area to collect the same ensquared 
energy given the Sloan telescope feed. 

16 After ~ 6 months of operation, 7 individual fibers within IFUs 
have broken (1 in manufacturing, 1 in assembly, 5 in operation), 
representing < 0.1% of the total. Detailed statistics on the break¬ 
age frequency of stripped, fused fiber bundles are unknown but 
would have represented a significant cost increase in manufactur¬ 
ing. 
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Fig. 2.— Effective IFU filling factor (live fiber core area divided 
by total IFU footprint) as a function of buffer thickness for an 
ideal 127-fiber (solid line) and a 19-fiber (dotted line) hexagonal 
IFU. The small difference between the solid and dotted lines repre¬ 
sents the diminishing importance of edge effects in the hexagonal 
footprint as the IFU area increases. The filled star represents the 
measured 56% filling factor of the as-built 127-fiber MaNGA IFUs 
IIDrorv et al.ll2015l ). which is consistent with theoretical expecta¬ 
tions. The filled triangle shows the 75% filling factor of the SAMI 
survey bundles (5 (im cladding) for comparison. 

Sydney-style bundles (whose fibers are chemically etched 
to an outer diameter of ~ 132 /im) achieve a filling fac¬ 
tor of ~ 70%, while the as-built hexagonal bundles with 
fully-buffered fibers achieve a filling factor of 56%. 

3. BASIC CONSIDERATIONS 

3.1. Required Performance 

Since the fiber bundles consist of 2” diameter circu¬ 
lar apertures separated by large interstitial gaps, each 
exposure will significantly undersample the point spread 
function (PSF) at the focal plane of the telescope (typi¬ 
cally ~ 1.5”) and produce a non-uniform response func¬ 
tion across the face of each IFU. We require that the 
MaNGA IFUs deliver sufficiently uniform performance 
that physical structures do not vary in shape as a func¬ 
tion of where they happen to fall within the IFU (i.e., a 
circular star forming region within a galaxy should ap¬ 
pear circular in the final MaNGA data cube regardless of 
whether it is in the center or the outskirts of the galaxy). 

A convenient way to place a limit on the level of unifor¬ 
mity required is to ensure that variations in the 2d PSF 
of the reconstructed MaNGA data cubes do not signifi¬ 
cantly impact measurements of the Balmer decrement or 
BPT-style (e.g JBaldwin et ahlll98lD line ratio diagrams. 
Since atmospheric differential refraction shifts the effec¬ 
tive position of each fiber as a function of wavelength 
Oi [O II] and Ha observations of a given Hu region 
for instance will be obtained with a slightly different con¬ 
figuration of fibers - while H a emission may be centered 
in a given fiber, spatially coincident [O II] emission may 
be centered in the interstitial region between fibers. 

As outlined in the MaNGA Science Requirements Doc¬ 
ument (SRD; see Yan et al. in prep), relative spectopho- 
tometry between [O II] (A = 3727 A) and Ha (A = 6564 
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Fig. 3. — Fractional error in the recovered flux from a point source 
if the assumed FWHM and axis ratio were incorrect. The refer¬ 
ence source is taken to have a circular gaussian PSF with FWHM 
2.5 arcsec; integrating the flux within a 2.5 RMS width aperture 
nominally encloses 95.7% of the total flux. If the actual FWHM 
is smaller (larger) than the model along any dimension the total 
flux enclosed by the aperture increases (decreases), resulting in an 
overestimate (underestimate) of the total flux. The dashed red 
line indicates the 7% error threshhold set by the MaNGA SRD; 
the solid black star indicates our adopted limits on the allowable 
variability of the delivered MaNGA PSF (15% in axis ratio, and 
15% in circularly averaged PSF FWHM). 

A) must be accurate to 7% or better in order to ob¬ 
tain the desired constraints on the star formation rate 
(SFR) and nebular metallicity within galaxies. We there¬ 
fore explore how this required spectrophotometric accu¬ 
racy translates to limits on the spatial variability of the 
MaNGA PSF. 

We begin by assuming that the PSF in a typical 
MaNGA reconstructed data cube can be characterized 
by a circular gaussian with a FWHM of 2.5 arcsec (as we 
discuss at greater length in Law et al., in prep, this model 
is a good approximation to the MaNGA commissioning 
data). Using typical aperture photometry techniques, a 
circular aperture of radius 2.66 arcsec (i.e., 2.5 times the 
radial scalelength of the PSF) would nominally enclose 
96% of the total fluxF 7 ! In Figure[5]we illustrate how de¬ 
viations from the nominal PSF model would affect this 
total; as the PSF becomes broader or more elongated 
the flux contained within the fixed aperture decreases, 
meaning that the derived aperture-corrected total fluxes 
would be in error{3 In particular, we find that an er¬ 
ror of 20% in the profile FWHM and 15% in the profile 
minor/major axis ratio is sufficient to bias the resulting 
flux measurements at the 7% level (filled star in Figure 
0 ). 

Similarly, in order to ensure that our limiting fluxes 
for undetected nebular transition features are accurate at 
the 7% level we also require that the signal-to-noise ratio 

17 If we were to adopt a PSF model with more power in the 
wings, or shrink the size of the circular aperture the variability 
between different PSF shapes would increase and lead to more 
stringent constraints on the allowable variability in the delivery 
MaNGA PSF. 

18 If the goal were to measure the flux from a single bright source 
whose structure is known a-priori to be effectively a point source 
then the actual light profile could be measured at each wavelength 
and the aperture adjusted accordingly. However, such a-priori 
knowledge of the intrinsic source structure cannot generally be 
assumed. Similarly, we assume that wavelength-dependant vari¬ 
ations from the A -1 / 5 Kolmogorov atmospheric turbulence profile 
are taken into account in determining the appropriate aperture. 










MaNGA Observing Strategy 


5 


of our data is constant at the 7% level across each IFU. 
Since the limiting flux is proportional to the square root 
of the exposure time, this translates to a requirement 
that the exposure time is effectively constant across each 
IFU at the 15% level. 

These three metrics (circularity, FWHM, and signal- 
to-noise ratio) therefore set our requirements on the uni¬ 
formity of the reconstructed image profile such that the 
calibrated fluxes derived from MaNGA data cubes are 
accurate at the 7% level. Ideally, however, we would 
prefer that spatial sampling issues not dominate the flux 
calibration accuracy budget for the MaNGA data cubes, 
and we therefore set a goal of achieving photometric per¬ 
formance at the 3.5% level where possible. The MaNGA 
hardware construction, dithering pattern, and observing 
strategy is therefore set by the following three high-level 
considerations: 

1. The reconstructed FWHM of all angular resolution 
elements in a bundle should vary by < 10% (goal) 
or 20% (requirement) across each IFU. 

2. The reconstructed minor-to-major axis ratio of all 
resolution elements in a bundle should be b/a > 
0.93 (goal) or 0.85 (requirement) across each IFU. 

3. The effective integration time of all resolution ele¬ 
ments in a bundle should vary by < 7% (goal) or 
15% (requirement) across each IFU. 

3.2. Integration Time 

The total integration time is set by our requirement 
that MaNGA reach a signal-to-noise ratio of 5 A -1 
fiber -1 in the r-band continuum at a surface brightness 
of 23 AB arcsec -2 . As described by Wake et al. (in 
prep.) and Yan et al. (in prep.) the typical integration 
time per plate to reach this target is anticipated to be 
about 3 hours in median conditions. In good conditions 
however the required time could be as low as 1.5-2 hours, 
and for particularly low-latitude fields the required time 
could be as much as 4-5 hours. This substantial varia¬ 
tion in total exposure time requires an observing strategy 
flexible enough to accommodate it. 

The optimal integration time for individual exposures 
is constrained by the MaNGA hardware and typical back¬ 
ground sky spectrum at APO. One of the strengths of 
MaNGA is the high throughput of the BOSS spectro¬ 
graphs shortward of 4000 A, and we therefore integrate 
each exposure for long enough that the shot noise from 
the background sky spectrum and detector dark current 
exceeds the read noise. The total noise N as a function 
of wavelength is given by 

N(X) = V(fsW + fdm)t + niN? (3) 

where / s (A) is the background sky spectrum in units of 
e - minute -1 per spectral pixel, / d is the dark current in 
e“ pixel -1 minute -1 , N r is the read noise in e - pixel -1 , 
nr = 3 pixels is the spatial width of a spectrum on the 
detector (see discussion by Law et al. in prep), and t 
is the integration time of an exposure in minutes. Re¬ 
arranging Eqn. [3] we find the time f m i n required for the 
combined sky background and dark current to equal the 
read noise: 


^min(A) 


Til IV 2 

/s(A) + / d ni 


(4) 


We estimate / s (A) for a typical MaNGA dark-time ob¬ 
servation using commissioning data from all-sky plate 
7341 (i.e., a calibration plate for which all IFUs target 
regions of blank sky) observed on MJD 566930 Fol¬ 
lowing the data model outlined by Law et al. (in prep), 
we take the FLUX array of the reduced mgFrame Hie 
(in units of flatfielded e - per spectral pixel), multiply 
by the SUPERFLAT array to obtain spectra in raw e - 
per spectral pixel, and combine ~ 600 individual fiber 
spectra to construct an extremely high-precision model 
of the background sky. We take the detector read noise 
to be R n = 2.0 (2.8) e - pixel -1 , and the dark current to 
be 0.033 (0.066) e - pi xel -1 minute -1 for the blue (red) 
camera (see Table 4 of lSmee et al.ll20i 3lF1 

We plot t m i n as a function of wavelength in Figure [H 
and note that the sky background rapidly dominates over 
read noise at almost all wavelengths, especially in the 
vicinity of strong OH atmospheric emission lines in the 
near-IR. The upturn in t min shortward of 4000 A repre¬ 
sents the falloff in blue sensitivity of the detectors, but an 
integration time of 15 minutes per exposure ensures that 
observations are shot-noise dominated for all A > 3700 

A. 

Although an integration time of longer than 15 minutes 
would further decrease the contribution of read noise to 
the total error budget, such longer integrations are un¬ 
desirable because of the cosmic ray event rate recorded 
by the red-channel detectors. In practice, the maximum 
integration time of each exposure is also limited by dif¬ 
ferential atmospheric refraction considerations (see £)7.2D , 
and we therefore adopt a nominal time of 15 minutes per 
exposure. Each completed plate will therefore consist of 
~ 6 — 20 exposures in order to reach the target depth. 


3.3. Numerical Simulations 
3.3.1. Simulation Method 

In order to assess the relative performance of differ¬ 
ent IFU bundles and observing techniques we perform a 
series of numerical simulations designed to test the uni¬ 
formity of their response to unresolved point sources (for 
which spatial structure is most pronounced). Adopting a 
working box size of ~ 45 x 45 arcsec with simulated pix¬ 
els spaced every 0.1 arcsec we first compute the footprint 
of a given IFU; this defines a mask image for which each 
fiber in the IFU is associated with a given set of pixels in 
the telescope focal plane that its light-sensitive core sub¬ 
tends. We then create an input ‘image’ to be observed 
by the simulated MaNGA IFUs by convolving a delta 
function by a model of the PSF at the focal plane of the 
SDSS 2.5-m telescope. This focal-plane PSF is taken to 
be the sum of two Gaussian profiles with FWHM 8 and 
2 8 respectively (where 8 = 1.4 arcsec is the FWHM of 
the median atmospheric seeing profile divided by 1.05) 

19 MJD (Modified Julian Date) 56693 corresponds to February 
5, 2014. 

20 The dark current is typically < 2% of the dark-time sky back¬ 
ground signal. 
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Fig. 4.— Exposure time t m i n required for a typical MaNGA dark-time sky spectrum to be dominated by Poisson noise from the background 
sky plus detector dark current. The break around 6000 A represents the dichroic break between red and blue channels; in reality there is a 
~ 300 A overlap between these channels. Strong features longwards of ~ 8000 A are due to bright OH sky lines. The dashed line indicates 
the adopted 15 minute exposure time. 


and peak amplitude ratio of 9/10 This input image is 
convolved with the top-hat fiber mask to determine the 
total amount of light received by each fiber; although 
the present simulation considers only a single input im¬ 
age the technique is immediately generalizable to multi- 
wavelength input image slices. 

We reconstruct a two-dimensional image from the in¬ 
dividual fiber fluxes using a flux-conserving variant of 
Shepards method (inverse-distance weighting) sim i lar t o 
that used by the CALIFA survey (iSanchez et all I2012T) . 
As part of MaNGA design simulations we explored alter¬ 
native methods of image recons truction such as drizzling 
(e.g., as adopted by SAMI, see lSharn et alJ 1201 5( 1. thin- 
plate-spline fits, minimum curvature surface fits, and 
kriging. As discussed by Law et al. (in prep) the modified 
Shepard’s method yielded the best results, and here we 
adopt the same parameters (e.g., final spaxel scale of 0.5 
arcsec) as used by the MaNGA Data Reduction Pipeline 
(DRP) for genuine survey data. The reconstructed image 
is fit with a 2d Gaussian model to determine its FWHM 
and axial ratio; major axis rotation is left as a free pa¬ 
rameter. 

This exercise is repeated for delta functions located 
in each of the 0.1 arcsec grid squares that lie within the 
central 75% of the IFU fiber bundle footprint (i.e., ignor¬ 
ing edge effects from point sources located on the outer 
ring of an IFU), resulting in ~ 40,000 simulated points 
across a 127-fiber IFU bundle. In Figure[5] (top row) we 
plot the on-sky footprint of a hexagonal fiber array, along 
with the variations in effective exposure time (exposure 
time multiplied by the fraction of the total light that is 
collected by fibers rather than being lost to interstitial 
regions), FWHM, and minor-to-major axis ratio of the 
reconstructed point spread function (PSF) as a function 
of the location of the point source within the fiber bun¬ 
dle. As anticipated, we note that all three quantities vary 
substantially across a given IFU in a single exposure. 

21 Mathematically, this is equivalent to the linear sum of 9/13 
times the input image convolved with a Gaussian of FWHM 9 plus 
4/13 times the input image convolved with a Gaussian of FWHM 
29. This profile provides a reasonable approximation of the on-axis 
SDSS focal plane PSF, matching the inner parts of the profile well 
and accounting for most of the flux in the outer wings (J.E. Gunn, 
priv. comm.). 


We quantify these results by calculating the RMS of 
the distributions in effective exposure time and recon¬ 
structed PSF FWHM (relative to the median values as 
[X - XmedianJ/Xmedian), the 3cr width W 99 encompass- 
ing 99% of these values, and the 99% lower bound for 
the minor-to-major axis ratio. For a single exposure, the 
effective in teg ration time varies by W 99 = 30.7% around 
the median 22 ! value; unsurprisingly, the greatest fraction 
of the total light is recorded for objects that are centered 
in a fiber, while the least amount is recorded for objects 
in interstitial regions. Similarly, the reconstructed PSF 
FWHM varies by almost 80% (from ~ 2 to 4 arcsec) de¬ 
pending on where a source falls with respect to the fiber 
grid, and the minor-to-major axis ratio b/a of the recon¬ 
structed image varies from ~ 0.5 — 1.0 (99% of values 
b/a > 0.53). 

In practical terms, this means that an unresolved H n 
region observed with such an IFU for just a single expo¬ 
sure may appear to be compact and circularly symmet¬ 
ric if it lands directly in the middle of a fiber, elongated 
and skinny if it falls directly between two fibers, or large 
and triangular if it falls midway between three adjacent 
fibers 0 Allowing for the effects of chromatic differential 
refraction ( jj4Tj) , this means that a single such H n region 
may simultaneously be sampled by all three different such 
configurations at different wavelengths. 

3.3.2. Dithering 

The sampling irregularities from fiber-bundle IFUs 
with substantial interstitial light l osses are well known 
from previous IF U surveys (e.g., I Sanchez et all 120121 : 
ISharp et all 120151 ). and can be largely overcome by ob¬ 
taining dithered observations. The geometry of the 
hexagonal fiber arrangement readily lends itself to a fixed 
triangular three-point dithering scheme that effectively 
fills the interstitial regions as illustrated in Figure [6] Re¬ 
peating the simulations performed in 1 13.3.11 with such 
dithered observations, we find that the combined data 
from just three exposures is able to achieve remarkably 

22 The median effective exposure time is just the filling factor 
(0.56) times the actual exposure time. 

23 Strictly, a single exposure simply does not have the spatial 
sampling in these cases to discriminate (for instance) between an 
unresolved point source and an elongated source. 
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Fig. 5. — Simulations of point-source response as a function of location within an IFU for a single exposure (top row) and a dithered set 
of exposures (bottom row) using a theoretically perfect hexagonal fiber bundle. The left-most panels show the footprint of the IFU fibers 
on the sky, the second column of panels show the percentage variations about the median exposure time as a function of position within 
the bundle. The third column of panels shows the deviation from the median delivered PSF, and the right-hand column of panels shows 
the recovered minor/major axis ratio. For undithered observations the greatest effective depth is obtained for sources located in the middle 
of a fiber (as is the smallest and most circular reconstructed image of a point source), while point sources falling in interstitial regions 
between fibers have minor/major axis ratios as low as ~ 0.5 and FWHM nearly double that of sources centered within a fiber. Numbers in 
panels 2-4 indicate the RMS deviation between values (rx), the 3cr width encompassing 99% of all values (W 99 ), and the minor/major axis 
ratio above which 99% of point lie (( 6 /a) 99 ). 


uniform image quality at all locations within a single IFU 
(Figure [5l bottom row). In contrast to the undithered 
case, 3-point dithering delivers effective exposure time 
constant to within 0.3% RMS, FWHM of 2.69 + 0.01 arc- 
sec, and ellipticity < 0.04. This uniformity easily meets 
the MaNGA science requirements described in 113.11 

Logically, the 3-point dithering pattern could be ex¬ 
panded to a regular 9-point pattern that also provides 
uniform coverage of the interstitial gaps, but with a finer 
sampling of the image plane. Although simulations sug¬ 
gest that this could provide ~ 10% improvement in the 
delivered PSF FWHM, such gains were not realized on- 
sky in tests with the MaNGA prototype hardware. This 
lack of improvement with respect to theoretical calcula¬ 
tions is likely due to a confluence of numerous complicat¬ 
ing factors, including degradation of the nominal dither¬ 
ing pattern by atmospheric refraction (see JH) , variations 
in fiber-to-fiber sensitivity, and changes in seeing and 
transparency conditions between exposures (see H6.2I) . 

3.3.3. Fiber Packing Regularity 

The gains achievable with such dithering depend fun¬ 
damentally on the uniformity of each IFU fiber bundle 
so that a single telescope offset can simultaneously dither 
each of our 29 IFUs (17 science and 12 calibration bun¬ 
dles) across the 3° field such that their fibers align with 
the interstitial gaps from the previous exposure. If fibers 
are not located at regular positions within every IFU, 
the dithering will not be able to uniformly sample the 
image plane. We explore the effect of fiber packing ir¬ 


regularity by repeating our earlier simulations with the 
introduction of a random perturbation to the position of 
each fiber in the simulated IFU fiber bundle, such that 
each fiber is slightly offset from its nominal position by 
some distance drawn randomly from a Gaussian distri¬ 
bution with a given rms. Each simulated IFU bundle is 
observed with a nominal 3-pt dither pattern as defined by 
Figure [6] Additionally, we simulate the effect of observ¬ 
ing the circular Sydney-style fiber bundle with a 7-point 
dither pattern (based on that adopted by the SAMI sur¬ 
vey) that compensates for the irregular fiber placement 
with greater filling factor and a larger number of dithered 
sampling points. 

As indicated by Figure [71 neither the dithered Sydney- 
style circular fiber bundle nor the 20 /im tolerance hexag¬ 
onal fiber bundles meet our target regularity goals, with 
a recovered PSF KW1 [JYp~1 varying by > 20% over the 
extent of an IFU (i.e., 2.66 + 0.12 arcsec with 99% values 
ranging from ~ 2.3 - 2.9 arcsec), and minor/major axis 
ratios as low as b/a ~ 0.8. In contrast, using a hexag¬ 
onal fiber array constructed to a tolerance of 5 /im rms 
with a 3-point dither pattern we expect to achieve a PSF 
FWHM that varies by less than 10 % ove r a given IFU. 

As detailed by iDrorv et al.l (120151) . the as-built 
MaNGA fiber bundles meet and exceed our target thresh¬ 
old with a typical fiber placement accuracy of 3 /im rms. 
Using the as-measured fiber metrologjo for 127-fiber 

24 We quote the average of the minor- and major-axis FWHM 
values. 

25 The final placement of individual fibers within an IFU can be 
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Fig. 6 .— Schematic diagram of the 7 central fibers within a hexagonally-packed MaNGA IFU, showing the 120 micron diameter fiber 
core and surrounding cladding plus buffer. The triangular figure shows the relative positions of the three dither positions; the fiber bundle 
is located at position ‘S’. The central (C) ‘home’ position is labeled, along with the north (N), south (S), and east (E) dither positions. 
The nominal plate scale of the SDSS telescope is 217.7358 mm/degree, or 60.48 microns/arcsec. 


MaNGA bundle ma024, we simulate the anticipated per¬ 
formance using this fiber bundle in row E of Figure [ 7 ] 
With a nominal dither pattern we expect to achieve a 
PSF FWHM which varies by less than 7% over an IFU 
(i.e., 2.66 ± 0.01 arcsec with 99% values ranging from 
~ 2.6 to 2.7 arcsec), and has a nearly circular profile 
everywhere with 6/a > 0.95. 

3.4. Baseline Observing Strategy 

The dithered observing simulations presented in 33.3.21 
and exposure time requirements described in 33.21 mo¬ 
tivate a nominal observing scheme in which targets are 
observed in sets of 3 dithered exposures (N-S-E) of 15 
minutes each. Given the regularity of fiber placement 
with each IFU and the locator pins that constrain each 
IFU to have the same position angle, correctly-dithered 
exposures can be simultaneously obtained for all IFUs on 
a given plate by simply offsetting the telescope pointing 
with respect to the guide stars. Since the coverage and 
image quality of a single set of three dithered exposures 
is known to be acceptably uniform, the total summed 
coverage of N such sets will also be uniform and have 
a depth of 0.75IV hours, allowing us to simply observe 
additional sets of 3 exposures until the combined data 
reaches our target signal-to-noise ratio of 5 A -1 fiber -1 
in the r-band continuum at a surface brightness of 23 AB 
arcsec -2 . 

Such a scheme provides us with considerable flexibil¬ 
ity to adjust our total exposure time in 45-minute in¬ 
crements without adversely impacting the delivered data 
quality whether there are 6 or 20 total exposures for a 
given galaxy. It is this flexibility as much as the dithered 
performance simulations themselves that drives us to 
adopt the regular hexagonal fiber arrays for MaNGA 

measured to an accuracy of better than 1 jim IIDrorv et alJ[2015 t 


rather than the SAMI-style circular fiber bundles, which 
rely upon a large number of exposures at many differ¬ 
ent dither positions to statistically fill in the interstitial 
gaps0 However, since this technique relies upon tightly 
controlling the fiber locations to provide uniform cover¬ 
age we must properly mitigate a variety of effects that 
will act to degrade this uniformity, and this goal in turn 
drives many aspects of the survey operation. 

4. ATMOSPHERIC REFRACTION 

As a photon passes through the Earth’s atmosphere it 
is refracted by variations in the density of the air. Un¬ 
der the usual assumption of a plane-parallel atmosphere 
with a vertical density gradient this bends the light from 
an astronomical target along the parallactic angle (the 
great circle connecting the target and the observers local 
zenith), causing astronomical objects to appear slightly 
higher in the sky than they truly are. Atmospheric re¬ 
fraction introduces significant optical distortions that ad¬ 
versely affect our ability to dither our IFU observations 
to the desired accuracy. Loosely speaking, the effects can 
be split into chromatic differential refraction and field 
differential refraction which we detail below. 

4.1. Chromatic Differential Refraction 

Atmospheric refraction is a function of atmospheric 
conditions (temperature, pressure, and relative humid¬ 
ity), zenith distance (i.e., the amount of atmosphere 
that an incoming photon must traverse), and wavelength. 
The impact of such refraction on astronomical observa¬ 
tions has been studied at some length i n the literature 
(e.g., iFilippenkolIl982t iCubv et aLlll998l and references 

26 Additionally, the hexagonal tapered ferrule construction tech¬ 
nique can be scaled up to bundles with large numbers of hexagonal 
‘rings’ without significantly degrading the packing regularity. 
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Fig. 7. — As Figure[5] but showing simulated point-source response as a function of location in an IFU for dithered observations of fiber 
bundles built to a variety of specifications. Row A simulates a Sydney-style 61-fiber bundle using a 7-point dither pattern. Rows B-F 
simulate a 3-point dither pattern applied to a hexagonal arrangement of 127 fibers with varying RMS deviations of each fiber from the 
nominal position (cr = 0 — 20 ^im). Note that for display purposes panel A is zoomed in slightly compared to panels B-F. In order to meet 
our uniformity criteria we require a < 5 ^m, which our as-built IFUs achieve (row E). 
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Fig. 8 .— Differential atmospheric refraction in arcsec of altitude 
relative to 5500 A for the MaNGA wavelength range as a func¬ 
tion of zenith distance. Calculations assume median conditions for 
APO with air temperature 10.5° C, 24.5% relative humidity, and 
atmospheric pressure of 730 mbar. 


Since differential refraction (particularly shortward of 
4000 A) can be substantial compared to the fiber radius 
of 1 arcsec (see Figure |SJ) the spectrum recorded by a 
single fiber is not strictly the spectrum of a single region 
in a given galaxy; it is a bent ’tube’ that traces differ¬ 
ent regions of the galaxy at different wavelengths. Most 
immediately, this means that the effective on-sky foot¬ 
print of the MaNGA IFUs can be shifted by up to ~ 1 
arcsec between blue and red wavelengths, requiring that 
the MaNGA data reduction pipeline (DRP; see Law et 
al. in prep) rectify the spectra to a common astrometric 
grid when reconstructing the data cubes. More prob¬ 
lematically, since the three exposures in a given dither 
set will be obtained at different hour angles the relative 
offset at a given wavelength will change between these 
three exposures and degrade the intended dither pattern 
coverage. 

At the guide wavelength of 5500 A, the three dithers 
will be executed properly. As illustrated in Figure [9] how¬ 
ever, at other wavelengths there will be variable shifts of 
the effective dithering pattern. These shifts can in some 
cases be comparable to the dither distances themselves, 
thereby degrading the effective dither pattern such that 
entire dither postions can be effectively ‘lost’ at certain 
wavelengths. As suggested by Figures [5] and [7] this pro¬ 
duces substantial and undesirable non-uniformities in the 
reconstructed image depth and recovered FWHM profile 
across the face of each IFU. 

4.2. Field Differential Refraction 


therein); we adopt estimates of the magnitude of refrac¬ 
tion r at a given wavelength relative to a fixed ‘guide’ 
wavelength developed by Enrico Marchetti for ESQ FI 
The direction of the refraction is along the local altitude 
vector for a given star; this corresponds to the parallactic 
angle 77 defined by the spherical triangle with vertices at 
the star, the celestial pole, and the local zenith. 


sin h cos cf> cos 8 

tan 77 = ———- — - (5) 

sin <p — sin o cos z 


where h is the hour angle (h > 0 towards the West), 
<f> is the local latitude (</> = 32°46'49" for APO), 6 is 
the target declination, 0 is the zenith distance cos z = 
sin <f> sin 8 + cos <j> cos 8 cos h, and 77 is defined in the range 


-180° to +180°. 

The SDSS 2.5-m telescope is equipped with an alt- 
az mount and all plates are observed with a position 
angle of 0 ° s o the amount of refraction in focal-plane 
coordinated is given by 


Axfocal = -r sin 77 ( 6 ) 

Aj/focai = -r cos 77 (7) 

i.e., at transit h = 0 , 77 = 0 °, and hence the entirety of 
the apparent refraction is along the yfocal direction Fn 


27 See http://www.eso.org/gen-fac/pubs/astclim/lasilla/difTrefr.html 

28 SDSS xfocal/yfocal coordinates are defined such that -fxfocal 
corresponds to -fright ascension and -fyfocal corresponds to 
-fdeclination. 

29 In the present work we neglect the relatively small effect of 
distortions introduced by the SDSS 2.5m optical system; these are, 
however, accounted for in the actual data pipeline described by 
Law et al. (in prep). 


In addition to varying with wavelength, both the mag¬ 
nitude and the direction of atmospheric refraction vary 
according to the location of an object on the sky, and 
the 3° SDSS plugplate field over which our IFUs are dis¬ 
tributed is sufficiently large that this variation cannot be 
neglected. As a given field rises, transits, and sets, the 
apparent locations of astronomical targets in the tele¬ 
scope focal plane shift. As described in SJ21 the SDSS 
telescope guider system compensates for this using guide 
fibers placed on astrometric standard stars distributed 
throughout a given field, and adjusts the overall shift, 
rotation, and scale of the focal plane to compensate. 
However, since field compression occurs along only a sin¬ 
gle direction (altitude) it cannot be fully corrected by a 
global change in the focal plane scale, leaving a resid¬ 
ual quadrupole term in the guider-corrected focal plane 
locations of the target galaxies (see Figure fl0l)F n l 

Such field differential effects are most noticeable when 
observing wit h single fibers or an array of slits (see, e.g., 
discussion bv lGubv et allllOfFi . for the 16’ x 16’ VIMOS 
field of view) since targets can rapidly shift out of the 
aperture. Hence, previous generations of SDSS that have 
used single fiber spectroscopy have been careful to ob¬ 
serve at hour angles close to that for which a given plate 
is drilled. In contrast, MaNGA is relatively insensitive 
to shifts in the effective centroid of an IFU since such 
shifts are small compared to the total field of view of 
each IFU (~ 30 arcsec for the 127-fiber IFUs)E3 The 

30 Field differential refraction is calculated using the SDSS plate 
design code located on the collaboration SVN repository. 

31 This effect is more important for the spectrophotometric mini¬ 
bundles which have a diameter of only 7.5 arcsec; as discussed by 
Yan et al. (in prep), large offsets of the spectrophotometric stan- 
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Coverage pattern: 5500 A Coverage pattern: 3500 A 



Fig. 9.— Illustrative figure showing degradation of the intended dither pattern due to chromatic differential refraction. In this example 
we assume a target at S = +60° was observed with a standard N-S-E dither pattern, but the three exposures were taken at hour angles 
h = — 4h, Oh, and +4h respectively (corresponding to parallactic angles 77 = —97°, 180°, +97°). The image on the left shows the offset due 
to chromatic refraction at 3500 A relative to the nominal center of a given fiber and defines the regularity statistic f 2 . While the achieved 
dither pattern is nominal at the guide wavelength (central panel), at 3500 A the fibers in positions S and E lie almost atop each other 
(right panel). 



Fig. 10.— Magnified illustration of the effects of field differen¬ 
tial refraction at the guide wavelength (~ 5500 A) across the 3° 
diameter SDSS plugplate. Black £ -b’ symbols indicate the nominal 
positions in focal plane coordinates of a randomly-selected set of 30 
target galaxies. These locations are computed assuming that the 
plate center has declination +7° and is observed at transit (h = 0 
hours); these correspond to the locations of the physically drilled 
holes in the plugplate into which the MaNGA IFUs are inserted. 
If the same plate were observed 4 hours later (h = +4 hours) the 
apparent locations of the galaxies in the focal plane would be dif¬ 
ferent due to field differential refraction both before (red asterisks) 
and after (green open boxes) guider corrections have been applied. 
Note that all offsets from the nominal positions have been magni¬ 
fied by a factor of 300 to enhance visibility; the maximum actual 
shift after guider corrections in this example is ~ 2 arcsec. 

MaNGA plugplates are therefore all drilled for transit 
(h = 0 hours), so the holes into which the MaNGA IFUs 
are inserted correspond to the expected focal plane loca¬ 
tions of the galaxies at this point in time. 

More important for MaNGA is the change in held dif¬ 
ferential refraction between exposures in a given dither 
set, which leads to degradation of the effective dither 
pattern akin to what was seen for chromatic differential 
refraction in Figure [9] As illustrated by Figure [10] the 

dard stars from the center of the calibration minibundles due to 
a combination of differential refraction, dither offsets, and other 
effects can complicate flux calibration. 


magnitude of this effect depends on the held declination, 
the hour angle h of exposures within a given set, and 
the location of an IFU within the plugplate. In the ex¬ 
treme example shown in Figure OH] (low declination, with 
exposures obtained many hours apart) the shift can be 
comparable to a hber diameter. In more realistic and 
typical cases (held center at 6 = +40°, observed at h = 0 
and h = +1 hours) the shift after guider corrections is 
typically < 0.1 arcsec. 

5. THE UNIFORMITY STATISTIC Cl 

Given the presence of both chromatic and held differ¬ 
ential refraction, no two exposures taken by MaNGA will 
have an identical hber sampling pattern even in the ab¬ 
sence of dithering. The primary driver of the MaNGA 
observing strategy is therefore mitigation of the impact 
of atmospheric differential refraction on the regularity of 
the dither pattern in order to achieve maximally-uniform 
data quality and depth within a given IFU. 

Given any two exposures separated by a time At there 
are vectors rl and rl defining the effective offset of a hber 
from its intended location on the target galaxy due to 
chromatic differential refraction, and si and si the offset 
due to uncorrectable held differential refraction effects. 
In our rectilinear focal-plane coordinate system the total 
shifts from differential refraction are given by: 


Aaq = ri sin 771 + x ■ si ( 8 ) 

Aj/i = rr cos rh+y ■ s[ (9) 

A ^2 = r 2 sin 772 + x • si ( 10 ) 

Aj /2 = r 2 cos 772 + y ■ si (11) 


where 771 and 772 are the respective parallactic angles for 
the two exposures and the vectors si and si are each 
projected into their components along the x/y focal plane 
coordinate system. The quantity of interest for survey 
planning purposes is the total distance between these 
shifted locations in the focal plane: 


U = V(Axi - Ax 2 ) 2 + (Ayi - Aj/ 2) 2 (12) 

In practice, we calculate U between the hrst and last 
exposures in a dithered set of three frames (see illus- 
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trative diagram in Figure mB Using the tools devel¬ 
oped in 1 13.31 we simulate four test cases where f2 ranges 
from 0 to 1”. We use the as-built MaNGA 127-fiber 
IFU ma024, and assume a standard three-point (N-S-E) 
dithering strategy in which exposure N is shifted by 12/2 
in the -Xfocal direction, exposure S is shifted by 12/3 in 
the -Yfocal direction and exposure E is shifted by 12/2 
in the +Xfocal direction (see, e.g., Figure [9]). Note that 
we are free to assume such symmetry because any shift 
common to all three exposures will simply result in a 
translation of the entire pattern. 

We show results for the expected exposure time, re¬ 
constructed PSF FWHM, and reconstruc ted ax is ratio 
uniformity as a function of 12 in Figure ms Figure 
□J suggests that so long as 12 < 0.2 arcsec observations 
should meet the target uniformity criteria outlined in 
1 )3. 3. 31 with FWHM 2.65±0.08 arcsec. At 12 = 0.4 arcsec, 
degradations in the reconstructed PSF uniformity and 
circularity start to become apparent; although the mean 
reconstructed PSF in the bundle has FWHM 2.65 ±0.14 
arcsec the total spread of FWHM values can be as large 
as ~ 0.3 arcsec, and 99% of locations have minor/major 
axis ratio greater than 0.85. By 12 = 1.0 arcsec the dither 
pattern is badly degraded, with reconstructed FWHM 
values varying by over an arcsecond depending on where 
a point source falls within the bundle. Our science re¬ 
quirements 1 113.II) therefore translate to a requirement 
that 12 < 0.4 arcsec, with the goal of reaching 12 < 0.2 
arcsec for the majority of observations so that it does not 
dominate the flux calibration accuracy budget. 

6 . MANGA OBSERVING STRATEGY 
6 .1. Set Lengths and Visibility Windows 

As described above, O is a complicated function of 
wavelength, integration time, target declination, hour 
angle, and location of an IFU on a given plate. How¬ 
ever, it is possible to define a series of relatively simple 
observing guidelines that will ensure that 12 stays below 
our 0.4 arcsec threshold. 

First, we note that 12 behaves nearly linearly with the 
amount of elapsed time between exposures in a given set, 
meaning that it is desirable to obtain all three exposures 
in the set as close in time to each other as possible. Since 
each exposure is 15 minutes long, we therefore require 
that all three exposures be obtained in a set length of 
one hour (i.e., the change in hour angle between the start 
of the first exposure and the end of the last exposure 
should be 1 hour or less, corresponding to 45 minutes 
between the effective midpoint of the first and last expo¬ 
sure). While we expect that each set of three exposures 
will typically last 48 minutes accounting for typical read¬ 
out times and overheads, this hour-long block provides 
necessary flexibility in scheduling, especially during vari¬ 
able weather conditions. 

We next calculate the expected 12 within a 1-hour long 
set as a function of the midpoint hour angle h se t of the 
set (h se t denotes the absolute value of the hour angle 

32 Each exposure is 15 minutes in length; we adopt the mid¬ 
point of each exposure as th e cha racteristic instant for purposes of 
calculating Q (although see <17.21 . 

33 Note that while £7 degrades the expect coverage pattern, we 
assume that the magnitude and direction of all of these shifts are 
known (see discussion by Law et al. in prep) and the true effective 
locations of each fiber are used when reconstructing the data cube. 


midway between the start of the first and end of the last 
exposure). In Figure fT2l we show the results of this cal¬ 
culation for three different wavelengths, three locations 
on a plate, and a range of different declinations[3 As 
expected, 12 is largest at extremely blue wavelengths (for 
which chromatic differential refraction is greatest) and 
on the edges of a plate (where uncorrected field differ¬ 
ential refraction is greatest). More importantly however, 
we note that 12 grows rapidly with increasing hour an¬ 
gle (either East or West of the meridian) meaning that 
we want to obtain our observations as close to transit as 
possible. Our 12 limit therefore equates to defining a se¬ 
ries of visibility windows around transit within which 
all MaNGA observations must be taken. 

In order to compute the length of these visibility win¬ 
dows we require that 12 must be less than 0.4 arcsec for 
all sets, at all wavelengths, at all locations on a given 
plate, and at all declinations. As indicated by Figure 
m the worst wavelength for 12 will be 3600 A, where 
the chromatic refraction is greatest. We work out the 
worst location on a given plate as a function of declina¬ 
tion by using Monte Carlo techniques to compute 12 for 
each of 20,000 randomly chosen locations on an SDSS 
plugplate over the course of a 1-hour set. As illustrated 
by Figure [13] the worst 12 is typically for IFUs located 
on the Eastern/Western edges of the plate for target 
declinations ~ +30 — 40°; this pattern shifts at more 
northerly/southerly declinations. 

Using these simulations we finally have all of the pieces 
required to define our visibility windows. For a grid 
of declinations spaced every 5° from 5 = 0° to 70° we 
compute the limiting set hour angle such that 12 = 0.4 
arcsec at A = 3600 A at the worst location on a given 
plate. Converting the set midpoint hour angle to the 
maximum midpoint hour angle of an individual exposure 
(2i ex p = h se t + 22.5 minutes for 1-hour sets), we show the 
final visibility windows as a function of declination in 
Figure 1141 These windows range from about lh either 
side of transit for fields near the celestial equator to ~ 3 
hours for declinations 6 ~ +40°. 

Intriguingly, despite all of the complications involved 
in computing these visibility windows they are nearly 
equivalent to simple airmass limits, independent of field 
declination. As illustrated by Figure [TT] our visibility 
windows can be described as a 6th order polynomial as 
a function of declination, or more simply by the require¬ 
ment that airmass AM < 1.21 for all exposures at all de¬ 
clinations. This airmass limit is determined by the SDSS 
plate diameter, the BOSS spectrograph wavelength cov¬ 
erage, and the assumed length of each set P’1 

We note that while these visibility windows have been 

34 Due to symmetries inherent in this exercise (chromatic and 
field differential effects combining constructively or destructively), 
at fixed wavelength one side of the plate will exhibit the worst at 
positive hour angles (west of meridian) and the other at negative 
hour angles (east of meridian). For convenience we collapse the 
problem such that h se t refers to the absolute value of the hour 
angle, and Q is taken to be the greater of the value from =bh S et- 

35 It is therefore possible to increase the airmass limit by reduc¬ 
ing the set length or effective plate diameter (i.e., restricting the 
locations of IFUs on the plate). For instance, a set length of 48 
minutes instead of 1 hour would increase the airmass limit to 1.34, 
expanding the visibility windows significantly. Such modifications 
to the observing strategy set forth here will be actively explored 
over the lifetime of the survey. 
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Fig. 11.— As Figure[5] but showing simulated point-source response variability as a function of location in an IFU for dithered observations 
of MaNGA 127-fiber bundle ma024 with different values of the pattern degradation Q. 


established to ensure that £7 < 0.4 arcsec at all wave¬ 
lengths for all MaNGA observations, typical performance 
is expected to be considerably better than this. At most 
wavelengths, most locations on a plate, and most hour 
angles within the visibility window £7 will be 0.2 arc- 
sec or below (see, e.g., Figs. H2l and H3l). Additionally, 
these simulations have assumed that sets are completed 
in one hour (45 minutes between the midpoint of first 
and last exposures in a set). Early survey observations 
at APO suggest efficiency such that most sets are actu¬ 
ally observed in more like 48 minutes (33 minutes be¬ 
tween the midpoint of first and last exposures); since £7 
scales roughly linearly with the set length we therefore 
expect on-sky performance to typically be a factor ~ 33% 
better than assumed in these simulations. Additionally, 
irregular coverage of an astronomical target in one set of 
exposures will tend to be averaged out across many such 


sets, resulting in more uniform performance for the final 
data cube of a given source. 

6 .2. Observing Conditions and Missing Exposures 

Thus far, all simulations have assumed that atmo¬ 
spheric seeing remains constant throughout all exposures 
in a given set, and that small variations in transparency 
can be normalized via per-exposure flux calibration (al¬ 
though see £17.41) . This assumption is often reasonable 
over the course of any given hour, but since rapid changes 
in observing conditions occur on some nights we must 
formulate our observing strategy accordingly. 

Consider, for instance, the pathological case where two 
dithered exposures have been successfully obtained in 
good conditions, but the third is lost. Whether it is 
never taken, or taken in extremely poor conditions (e.g., 
heavy cloud, seeing greater than 4 arcsec FWHM, etc.), 
the combined set of exposures no longer uniformly sam- 
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Fig. 12.— SI statistic as a function of midpoint hour angle of the set (h se t) for a range of wavelengths, target declinations, and locations 
on a plate. Left, middle, and right columns respectively show results for wavelengths of 3600, 5500, and 9000 A; top, middle, and bottom 
row respectively show results for an IFU in the middle of the plate, 1.06° towards the E edge of the plate (a circle at this radius encloses 
50% of the plate area), and on the E edge of the plate. Red, orange, green, blue, black, and grey solid lines respectively indicate results 
for declinations <5 = 0°, 15°, 30°, 45°, 60°, and 75°. The horizontal dotted lines at 11 = 0.2 and 0.4 indicate the thresholds of ideal and 
acceptable performance respectively. High-frequency structure in some lines is due to discrete changes in the best-fit guider corrections 
between individual simulation points. 


pies the source image. In such a situation, the missing 
exposure would have to be made up on another night, 
and obtained within a small range of allowable hour an¬ 
gles such that the total set length is still less than one 
hour. 

We therefore establish a series of additional require¬ 
ments for image uniformity across exposures within a 
given set. Based on simulations similar to those de¬ 
scribed in 1 13.31 and )J5] but allowing for variable seeing 
and transparency, we find that 

• All exposures in a set should have seeing within 0.8 
arcsec of each other. 

• All exposures in a set should have (S/N) 2 values 
within a factor of two of each other. 


• Each set of exposures should have median seeing 
2.0 arcsec or below in order for the reconstructed 
image to have FWHM less than 3 arcsec (ensuring 
uniformity of image quality between galaxies in the 
MaNGA survey). 

Historical conditions at APO and experience during 
MaNGA commissioning suggest that atmospheric con¬ 
ditions are generally stable enough that these criteria 
will not pose a serious limitation to survey operations. 
In practice, exposures also can often be rearranged be¬ 
tween sets to optimize observing efficiency and minimize 
the need for patching of missing dither positions (see 
discussion by Yan et al. in prep), and further modifica¬ 
tions to the baseline strategy will continue to be explored 
throughout the survey. 
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Fig. 13.— as a function of location on a plate centered at 
S = 40°. Simulations are performed at 3600 A and assume a 5- 
hour observing window (i.e., h S e t = 2.5 hours either side of transit). 
Each point represents the maximum value of D experienced at a 
given location for a hour-long set of exposures taken within this 
observing window (for one side of the plate this maximum will 
occur prior to transit, for the other side it will occur after transit). 



Declination (degrees) 

Fig. 14. — Black asterisks show the maximum hour angle away 
from transit (h exp ) within which all MaNGA exposures must be 
obtained as a function of declination based on numerical simula¬ 
tions. The solid black line represents a polynomial fit to these 15 
data points. Dotted lines indicate contours of constant airmass 
(every 0.05 from AM 1.05 to 1.30) as a function of declination and 
hour angle; note that these contours closely track the derived hour 
angle limits. 

7. ADDITIONAL CONSIDERATIONS 

Although differential refraction considerations are the 
primary factor that sets the MaNGA observing strategy, 
we also highlight a few additional considerations here 
that will impact the MaNGA reconstructed image quality 
and must be accounted for in survey operations. 

7.1. Required Dithering Accuracy 

Just as differential refraction effects degrade the ef¬ 
fective dithering pattern and contribute to non-uniform 
sampling of an astronomical source, so too does the 
dithering accuracy of the telescope. As described above 
in m U from refractive sources will frequently be less 
than 0.1 — 0.2 arcsec, and the individual telescope off¬ 
sets must therefore be good to better than 0.1 arcsec 


in order to not be the limiting factor governing the im¬ 
age sampling regularity for the majority of observations. 
Indeed, it is particularly important to minimize the con¬ 
tribution of offsetting errors for cases with already-high 
Q from differential refraction as the compounded errors 
may easily make the difference between an acceptably- 
versus unacceptably-uniform set of exposures. Based on 
observations performed at APO during MaNGA commis- 
sioningH the dither offset error has a median of 0.063 
arcsec, and is smaller than 0.1 arcsec in 76% of expo¬ 
sures. Although the current dithering accuracy degrades 
to a median of 0.1 arcsec at altitudes higher than 80°, 
work is ongoing to improve this performance (see details 
in Yan et al. in prep). 

7.2. Required Guiding Accuracy 

In addition to the accuracy with which the telescope 
offsets are performed it is also important to consider the 
guiding accuracy of the telescope (i.e., how well a given 
position is maintained over the course of an exposure). 
Although poor guiding performance will not degrade the 
coverage uniformity of a set of exposures, it will degrade 
the image quality of the exposures by contributing in 
quadrature to the effective astronomical seeing. Obser¬ 
vations obtained during MaNGA commissioning show 
that the median guiding accuracy (based on variations 
in guide star positions across all 15-second guider cam¬ 
era exposures during each 15-minute science exposure; 
see details in Yan et al. in prep) is 0.12 arcsec, substan¬ 
tially smaller than the median SDSS 2.5m seeing of ~ 
1"5 (computed across all BOSS spectroscopy in 2012). 

We note that a similar effect is caused by differential 
refraction; just as changing refraction causes the effective 
location of a fiber to move between two exposures (1 14. 21) . 
so too does it cause the effective location of a fiber on a 
given astronomical target to move during the exposure 
as well. However, for observations obtained using the 
strategy outlined in f[6] above this effect is small. Since 
differential motion of a fiber with respect to a fiducial 
position (i.e., D) scales roughly linearly with time, the 
motion in a 15-minute exposure will be ~ 1/3 of that for 
a given set of 3 exposures. Since we require the latter 
to be < 0.4 arcsec even in extreme cases, the motion in 
a single exposure will be < 0.1 arcsec, which is small 
compared to the guider accuracy and atmospheric seeing 
profile. 

7.3. Bundle Rotation 

The rotational position 8 of the MaNGA IFU bun¬ 
dles is controlled via clocking pins that ensure proper 
alignment of each IFU. However, mechanical tolerances 
of the pinhole translate to a rotational uncertainty for 
each bundle at the level of ~ 3°. In order to ensure that 
individual sets meet the U < 0.2 arcsec coverage regular¬ 
ity goal at the edges of the largest fiber bundles (~ 16.5 
arcsec radius) we specify that the rotational offset Ad 
required between any two exposures in a given set be 

A9 < tan _1 (0.2/16.5) = 0.7° (13) 

36 The guider system uses 16 coherent imaging fiber bundles 
plugged on the plate and imaged by a separate guider camera; by 
monitoring the positions of these 16 stars and comparing them 
to the desired positions we can measure the dithering accuracy 
directly. 
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Generally, rotational tension in the IFU cables should 
ensure that 9 remains relatively constant for a given plug¬ 
ging, and preliminary tests indicate that Ad ss 0.2° (see 
Law et al, in prep for further details). However, changes 
in the routing path of each IFU cable through the car¬ 
tridge can lead to a different rotation (and small trans¬ 
lational offsets) each time the plate is replugged, and we 
therefore require that sets be completed within a single 
plugging. 

Rotation between sets of exposures can also degrade 
the reconstructed image quality if it is not measured and 
accounted for in the final astrometric solution. In Fig¬ 
ure ED we simulate the effect of stacking together two 
sets with different rotations without accounting for their 
rotational offset in the fiber astrometry. Visible degrada¬ 
tion of the reconstructed PSF starts to become apparent 
at the edges of the IFU bundle once the sets are rotated 
from each other by ~ 3°, and distortions become severe 
once the rotation reaches 6° (i.e., ~ a fiber radius at the 
edge of the largest bundles). We therefore require that 
the MaNGA data pipeline be able to measure the rota¬ 
tional clocking of each IFU at the level of ~ ±1° so that 
it can be incorporated into the astrometric solution. 

7.4. Errors in Spectrophotometry 

As discussed by Law et al. (in prep) and Yan et al. 
(in prep), each MaNGA exposure is flux calibrated in¬ 
dependently to account for variations in the atmospheric 
seeing and transparency. Adequate image reconstruction 
is therefore dependent on the relative accuracy of the flux 
calibration between exposures in a given set; any offsets 
between exposures will hamper the ability of the dithered 
exposures to properly sample the source profile. The 
most pronounced effect of such offsets is not their degra¬ 
dation of the spatial profile of unresolved structures (e.g., 
point sources) however, but rather their introduction of 
artificial spatial structure into a smooth background. 

We therefore simulate dithered observations of a con¬ 
stant surface-brightness source (e.g., like the outskirts of 
a smooth elliptical galaxy), assuming typical observing 
conditions with visual seeing ~ 1.5 arcsec. We mimic 
flux calibration errors by multiplying the fiber fluxes for 
each exposure by a scale factor drawn randomly from a 
gaussian distribution with a given RMS width and a me¬ 
dian of 1.0 before reconstructing the composite image. 
As illustrated in Figure M calibration errors between 
individual exposures results in a stippling of the smooth 
background, introducing artificial spatial structure cor¬ 
related with the dithered fiber pattern. 

In a single set of 3 exposures, we find that RMS flux 
calibration errors of 2% between exposures results in a 
reconstructed image whose surface brightness varies by 
0.4% (RMS) from pixel to pixel. This is comparable to 
the 0.3% pixcl-to-pixel variations that we find in the re¬ 
constructed image assuming perfect flux calibration of 
all exposures. As flux calibration accuracy degrades fur¬ 
ther to 5%, 15%, and 50% RMS between exposures, we 
find pixel-to-pixel variations of 1%, 2%, and 8% respec¬ 
tively in the reconstructed image. If flux calibration er¬ 
rors are uncorrelated between exposures in different sets, 
this variation will average out over the course of obser¬ 
vations for a given plate. Even in the case where indi¬ 
vidual exposures are calibrated as poorly as to within 
a factor of 2, the pixel-to-pixcl flux for a uniform back¬ 


ground source varies by just 2% RMS when averaged 
over 4 sets (12 total exposures). In contrast, preliminary 
results from MaNGA commissioning data indicate that 
individual exposures are typically calibrated to within 
2.5% (Yan et al. in prep), suggesting that flux calibra¬ 
tion errors are unlikely to contribute significantly to the 
image reconstruction fidelity. 

8 . SUMMARY 

The MaNGA hardware design and observing strat¬ 
egy is driven by the desire to ensure high, uniform im¬ 
age quality and depth across all 10,000 of the galax¬ 
ies that will be observed during SDSS-IV. In particular, 
the goal of reaching 7% spectrophotometric accuracy be¬ 
tween [O II] A3727 and Ha A6564 requires that the recon¬ 
structed PSF varies by 10% or less (in both width and 
cllipticity) across the face of each IFU. This goal is par¬ 
ticularly challenging given the variable total number of 
exposures per field (~ 6 — 20) required to reach the target 
depth, chromatic differential refraction arising from the 
large wavelength coverage of the survey (AA3600—10,000 
A), and field differential refraction caused by the 3° wide 
field of view over which individual IFUs are deployed. 

We summarize the requirements necessary to meet our 
goal as follows: 

1. Each IFU fiber bundle should be constructed of a 
regular hexagonal grid of fibers to an accuracy of 
5 pm rms. The MaNGA IFUs meet and exceed 
this specification with a filling factor of 56% and a 
typical fiber placement accuracy of ~ 3 pm rms. 

2. Exposures should be obtained in sets of three 15- 
minute exposures dithered to the vertices of an 
1"44 equilateral triangle in order for each set to 
uniformly sample the image plane. 

3. The telescope must be able to dither to an accuracy 
of 0"1 or better. 

4. Each plate should be observed for an integer num¬ 
ber of sets until the combined depth reaches a 
signal-to-noise ratio of 5 A -1 fiber -1 in the r- 
band continuum at a surface brightness of 23 AB 
arcsec -2 . 

5. All three exposures in a set must be observed 
within one hour of each other (i.e., the change in 
hour angle between the start of the first exposure 
and the end of the last exposure should be one hour 
or less), and in a single plugging of a given plate. 

6 . All three exposures in a set should have (S/N) 2 
within a factor of two of each other, and be ob¬ 
tained in atmospheric seeing that varies by less 
than 0'/8. Each set should be obtained in median 
seeing of 2'/0 or better. 

7. All MaNGA exposures should be obtained in vis¬ 
ibility windows ~ 2 — 6 hours in length corre¬ 
sponding to airmass < 1.21. 

8 . MaNGA relative flux calibration between expo¬ 
sures must be good to ~ 5% or better. 


MaNGA Observing Strategy 


17 



0 0.5 1.0 -20% 0% +20%-10% 0% +10% 0.7 0.85 1.0 

Areal Footprint Relative Exp Time Relative PSF FWHM Axis Ratio 


Fig. 15.— As Figure [5] but showing the impact of uncorrected bundle rotation on the reconstructed image quality. Note the marked 
increase in FWHM and ellipticity of reconstructed point sources near the edges of the bundle for offsets ~ 6°. 


In reality, many of the issues considered here will tend 
to average out over the course of the 3-4 sets that will 
typically be obtained on a given plate since effects that 
cause a slight gap in coverage for one set will often be 
filled in by another. However, our objective in designing 
the MaNGA program is to ensure that the depth, cov¬ 
erage, and image quality of the entire survey is as good 
and uniform as possible for the entire wavelength range 
of each of our 10,000 galaxies. In forthcoming contribu¬ 
tions (Law et al. in prep, Yan et al. in prep) we will 
explore in greater detail how well we succeed in meeting 
these goals. 
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Fig. 16. — Simulated MaNGA observation of a constant surface- 
brightness field showing the characteristic stippling pattern intro¬ 
duced by relative flux calibration errors between dithered expo¬ 
sures. Such errors introduce artificial spatial structure correlated 
with the dithered fiber pattern. Greyscale stretch is arbitrary; 
in this example the blackpoint (whitepoint) is set to 20% below 
(above) the mean flux, corresponding to an RMS variation of about 
8% over the field of view. Note that the simulated field has been 
trimmed to omit effects from regions at the edge of the IFU field. 
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